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The Smoluchowski equation with a time dependent delta function sink is solved exactly for many 
special cases. In all other cases the problem can be reduced to an integral equation. It is shown that 
by knowing the probability distribution at the position of sink, one can derive analytical expression 
for probability distribution everywhere. Thus the problem is reduced from a PDE in two variables 
to an integral equation of one. As far as the authors knowledge, we are the first one to provide an 
exact analytical solution of Smoluchowski equation for a parabolic potential with time dependent 
sink. 


Diffusion-controlled reactions have been classical topics [if. There is huge interest in this subject because of it’s wider 
connections recognized in physics, chemistry and biology Q. One particular model of this class is a diffusion under 
a random-trapping environment. This problem can be modeled by a particle undergoing diffusive motion in presence 
of a time dependent Delta function sink. One may describe the diffusion motion by Smoluchowski equation and the 
sink could be a Dirac Delta function of arbitrary position and strength. As per the authors knowledge this Letter is 
the first one to consider the exact analytical solution of Smoluchowski equation with time dependent sink. All the 
solutions available in literature are the cases where the strength of the Dirac delta function used to model the sink 
is assumed to be time independent although the sink has to be time dependent in general in most of the cases 

and but the way strength of the sink depends on time is not well understood. In contrast, in the present work we 
have considered the cases where the strength of the sink varies with time. There are two common approaches one 
can think of using for solving this type of problems. One is based on using path integral method of Feynman type 
and the other is based on Laplace transforms. We use the latter method, although both are closely related. In the 
following we would like to solve the Smoluchowski equation with a time dependent sink. 

^P{x,t) = [L - k{t)d{x)]P{x,t). (1) 

In the above we denote that the probability of the particle to survive on the potential energy surface as P{x,t), and 

^ .9^ Ad dVix) 

The schematic view of our problem is given in Fig. I. We define the Laplace transform P{x,s) = L[k{t)P{x,t)] = 
P(a;, then the Laplace transform of Eq. (1) can be written as 

[s - L)] P{x, s) + L[k{t)S{x)P{x, t)] = Po{x, 0). (3) 

We will now solve Eq. (2), however we must first consider solving the following homogeneous equation Q 

[s - L]P{x,s) = Po{x,0). (4) 


Solution of this homogeneous equations must satisfy the boundary conditions given by 




In case of a parabolic potential represented by V(x) = Eq. (4) can be written as 
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P{x, s) = Po{x, 0) 


( 6 ) 


with P = Now we introduce the variable z and zq by z = x{^)^,Zo = xo{A/. The solution of 
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This, when substituted back into Eq. (9) gives 
P{z,s) = - 


ko DA0)D4mi - v)[^]^ fZ 


A 


1 + 


fco z>„(o)DUo)r(i-^)[^] 


i?:(-^<)i?J(-^>)exp[(zg - z^)/ 4 ]r(l-v)l^]i 


P(zo)dzo. 


Also in some situations it is always advantageous to present the solution in the following form 

/ OO 

G(z,sjzo)P(zo)dzo, 

-OO 

where G{z,s\zq) is nothing but Greens function. By comparing Eq. (10) and Eq. (11) we find 


( 12 ) 


(13) 


G{z,s\zo) = - 


ko i:i„(o)i:)„(o)r(i - u;(o)g;(o)exppgy4]r(i 


A 


1 + 


fcn g;(o)g;(o)r(i-G[?g 4 ]' 


DU-z^)D*{-z>) exp[(z2 _ z^)/4]r{l - v)[^]^ 


(14) 


Our expression for G(z, s|zo) as given above is the same as the one derived by Sebastian using another method 
Now we consider the case, where k{t) = —at, so that an equation similar to Eq. (8) would be 


dP{0+,s) dP{0-,s) 


dx 


dx 


T 

= -a{s) = -[k{t)P(0,t)] = I e-"' - ^P{Q,t)dt 


t 


( 15 ) 


Putting this solution for a(s) in Eq. (7) we will get 


P{z,s) = - 


A 


_5 

ds 


P(0,s) 


Dp(-z<)Op(-z>) exp[-zV4]r(l - v)[^] = 


Dl{-z<)Dl{-z>)eyi^[{zl - z^)/4]r(l - t;)[^] 2 


P{zo)dzo. 


(16) 


Now we consider P{zo) = S(zq — zq) as the initial condition, then above equation simplifies considerably to the 
following one (after we replace Zg by zq) 


a 


= -j 




Dy(-z<)Dy(-z>)expl-z^/4]r(l - v)[^] 2 


Dl{-z^)Dl{-z^)exp[{zl - z^)/4]r(l - v)Mi 


(17) 


For z = 0 the above equation can be rewritten as 


P{0,s) = -j 


D 40 )D 40 )r{i-v )[^]-2 


i4:(0)i^:(0)exp[(^2)/4]r(l-^)[^]i 


(18) 


Dividing the above equation on both sides by iM2HZ£Mp(i _ I ^ get 


D,( 0 )D,( 0 )r(l-u)[ 


2irAJ 

P(0,.) = -5^7(0,.) + g:<°)A)(0)exp|.g/4| 

27rA^ ^ ^ 


D„(0)D„(0) 


( 19 ) 
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Let us consider the following substitutions, we put 


/3(s) = 


D 40 )D 40 )r{l-v )[^]2 


and 


7(s) = 


i4;(0)j;(0)exp[^g/4] 

D,{0)D,{0) 


Hence Eq. (19) can be rewritten as 


ap,o,,) + M£)p(o,.) = 4 lW. 

as a a 


( 20 ) 

( 21 ) 

( 22 ) 


Multiplying the above equation on both sides by where and using some simple mathematics we 

reach at finally 


P{0,s) = 




7(s')ds', 


(23) 


^fg/3(s')ds'. _ 

able to determine the probability distribution function P{z,s) everywhere using Eq. (16). 


where /(s) = P{s')ds'. Hence we have found exact analytical solution for P(0, s) and once it is known we are 


Now we consider the case, where k{t) = —f, so that an equation similar to Eq. (8) would be 


dx 


dP(0+,s) dP{0-,s) . M /■“ . X, 

-dx-^ ^ = J e -P{0,t)dt 

/ OO pOO 

P{0,s')ds ; a{s) = jJ P{0,s')ds . 


(24) 


Putting this value of a{s) in Eq. (7) we get 


P{z,s) = - 


P(0, s')ds 


Dy{-z<)Dy{-z>) exp[-2V4]r(l - v)[^] 2 


D*{-z<)D*{-z>)exp[{z^ - z2)/4]r(l - v)[:^] 2 


P{zo)dzo. 


(25) 


Now we consider P{zo) = 6 {zq — zq ) as the initial condition, then above equation simplifies considerably to the 
following one (after we replace Zq by zq) 


P{z,s) = 


A 


P(0, s')ds 


Dy{-z<)Dy{-z>) exp[-zV4]r(l - v)[^] 2 


s 

D*v{-z<)Dti-z>) exp[{z^ -z^)/4]r(l -u)[^]^ 


(26) 


Eor z=0 the above equation can be rewritten as 


P(0,s) = 


A 


P(0,s')ds' 


p„(o)p„(o)r(i-u)[^]^ 


p;(o)p;(o)exp[zg/4]r(i -u)[^]^ 


Dividing the above equation on both sides by iM£^Ml!lr(i _ 2 ^ we get 


a 


p„(o)p.(o)r(i-u)[^]r H 


-P(0,s) = -/ P(0,s')ds + 


P:(0)P);(0)exp[zoV4] 


B 40 )Dy( 0 ) 


(27) 


(28) 


Hence the above equation can be rewritten as 

/ OO 

/ 

P(0,s')c^s +7(s)- 


(29) 
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In solving the above equation, we put u{s) = P{0,s')ds , then = P(0,s), hence above equation can be 

further written as 


du{s) a 7 (s) 

Ms) = 


ds M{s) 


PM' 


(30) 


Multiplying both side of the above equation by where f'{s) = — M{s) using some simple mathematics we 

reach at finally 


which can be used to evaluate P(0,s), as given below 


P(0,s) = 


ds 




PM) 


(31) 


(32) 


where /(s) = —^ fg /0J- Hence we have found exact analytical solution for P(0, s) and once it is known we are able 
to determine the probability distribution function P(z, s) everywhere using Eq. (26). 


Now we consider the case, where k(t) = Pexp{—at), so that an equation similar to Eq. ( 8 ) would be 

= -a(s) = 4 [k{t)P{Q, t)] = f e-"‘/3exp(-at)P(0, t)dt 

CLX J\. In 


dx 


Putting this value of a(s) in Eq. (7) we will get 




P(z .,) = -4p( 0,, + P.(-»<)P.(-»>)exp|-»V4]r(l - v)l^]z 

A s 


DvMMDvMz>)(iMi4 - 


P{zg)dZg. 


(33) 


(34) 


Now we consider P{zg) = 6 {zq — Zg) as the initial condition, then above equation simplifies considerably to the 
following one (after we replace Zg by zg) 


P{z, s) = s + a)DyMz<)DvMz>) exp[-2:^/4]r(l - ^ 

+D*vMz<)Dl{-z>) exp[(z2 _ z2)/4]r(l - ■ 


(35) 


For z=0 the above equation can be rewritten as 


P( 0 , s) = -4 p(0, s + a)P„(0)P„(0)r(l - v)[^]^ 


+P:(0)P:(0)exp[zoV4]r(l-u)[. 


B 


(36) 


^2TrA‘ ■ 

Further to obtain a solution of the above equation we could iterate the expression repeatedly to obtain the series 
solution, given by 


-f’(0:s) = ManM- 


(37) 


n=0 


Since we already have the form of this series, we may substitute Eq. (37) directly to Eq. (36) and then solving for 
a (s) by equating like powers of 77 , we find that 


P^(0)P^(0)exp[z^/4]r(l-u)[^]2 
“ As 


a„(s) = - 


P„(0)P„(0) exp[2g/4]r(l - v)[Ma] " 

■p*(o)p:(o)r(i-u)[^]^' 

Ans 

As 


for n = 0 
for n > 0 


(38) 
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It was shown in this letter that the Smoluchowski equation for parabolic potential with a time dependent Dirac 
delta function sink can be solved exactly in many special cases. Our method can be used to find exact solution to 
the case where the sink is represented by a non-local operator, and may be represented by \f > ko < g\, where / and 
g are arbitrary acceptable functions. Choosing both of them to be a Gaussian would be an improvement over the 
localized sink model. Sink may even be represented by a linear combination of such operators. 
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